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ABSTRACT 



SEA PLOT is a three-dimensional graphics program written 
for the Compass Hammer research effort. The program portrays 
a war-at-sea scenario in which a single aircraft attacks a 
single ship target located at the origin of an earth-fixed 
coordinate system. To plot the three-dimensional scenario 
on a two-dimensional plotting surface, a conic projection 
technique is used in which a cone of vision is generated 
based on the position of a viewer's eye. The position of 
the eye is automatically placed based upon the initial 
coordinate position of the attacking aircraft. The program 
is written using VERSATEC software available on the NPS 
IBM 3033 computer. However, the program has been structured 
in such a manner as to permit easy transition to a more 
sophisticated, picture oriented, graphics language such 
as PLOT 80 which would permit an interactive graphics 
employment of the program and possibly animation of the 
attacking scenario. 
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I. INTRODUCTION 



This thesis provides a program description and user's 
manual for the graphics package SEA PLOT designed for use 
in the Compass Hammer research effort. SEA PLOT is a 
three-dimensional graphics program which provides the user 
a perspective display of a war-at-sea scenario in which 
a single ship target is attacked by a single aircraft. 

The plot depicts (1) ship position and heading, (2) air- 
craft flight profile including coordinate position, roll, 
pitch and yaw, and (3) ship defensive missile launch envelope. 
All information is displayed in a three-dimensional format 
using a conic projection in which an optimal viewing 
position is automatically selected based upon the attacking 
aircraft's initial coordinate position. Because the plotting 
package is intended to provide a visual overview of the 
attacking scenario, it is not a reliable source for making 
measurements of ranges, bearings, etc. This results from 
the conic projection technique used which displays solid 
objects in three -space much as they would be viewed by the 
viewer's eye. Consequently, range marks and rings located 
in close proximity to the viewing position will show a 
wider dispersion than those at a more distant location. 

The conic projection itself is discussed in detail in 
Chapter II. 



Plotted output from SEA PLOT consists of range rings 
from the center of the target, aircraft flight profile, 
target ship orientation, and missile launch envelope. 

Range rings are displayed at 5000 -meter increments from 
the center of the target located at the plot center. 

Maximum range is calculated by examination of the attacking 
aircraft's flight profile such that maximum plot range 
corresponds to the maximum separation between target and 
attacker. The target ship is shown in planform only, but 
is oriented along a user input heading. Three-dimensional 
display of the target ship, which is possible with only 
minor program changes, resulted in a cluttered plot output 
and consequently was omitted. Display of a defensive 
missile launch envelope consists of a cylinder comprised 
of unconnected dots corresponding to maximum missile slant 
range. The cylinder is truncated at a height of twice the 
maximum aircraft flight profile altitude. Aircraft ground 
track is depicted via a series of dark dots plotted on the 
surface corresponding to aircraft position over ground. 

The air track is shown in perspective view as a bold line 
tracking the vehicle movement through the air. Along the 
air track, an aircraft body is positioned and drawn 
corresponding to a minimum of 2000 meters of ground travel. 
This body plot displays aircraft maneuvers in roll, pitch, 
and yaw. Aircraft altitude perspective is enhanced via 
vertical dashed lines connecting ground and air tracks. 
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SEA PLOT also contains a variety of input scaling 
parameters which permit the user to focus on the close-in 
encounter between the target and aircraft, or to expand 
the view and observe the entire profile from initial to 
final aircraft positions. This is accomplished by entry 
of the scale parameter, one of the program data input 
elements. For example, SCALE=1 causes a wide angle display 
of the scenario. SCALE=2 in essence magnifies the plot size 
by a factor of two, focusing on the area from the origin 
out to one half the original range. Figures 1-1, 1-2, and 
1-3 contain examples of plot output using SCALE=1, 2 and 
3, respectively. Chapter III discusses format and the 
required program inputs for SEA PLOT. 

Because SEA PLOT output is to be included in the Compass 
Hammer report, hard copy, printed output was required, 
as opposed to an interactive graphic output onto a CRT 
terminal. Consequently, the plotting package was written 
using VERSATEC software available on the NPS computer. 
Because this software system is oriented primarily toward 
display of data in graphic or tabular form, and not toward 
the drawing of pictures , the programming algorithms were 
more complicated than would be required for an interactive 
graphics software package, such as PLOT 80 used on many 
Tektronix machines. However, because the program does have 
many interactive graphics applications, such as an 
iterative design of an aircraft mission profile, it has 
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been structured as much as possible to permit a relatively 
simple conversion to other software systems. To accomplish 
this, the graphics instructions are written in internal 
subroutines, each of which contains the commands peculiar to 
the particular software system employed. For example, to 
draw a line connecting an array of coordinate points, 

SEA PLOT calls an internal subroutine DRAW. Resident in 
DRAW are the VERSATEC commands necessary to draw the line. 

To convert to another graphics system, one would merely 
have to adjust the instructions in DRAW to facilitate the 
software conversion rather than search through the entire 
program locating specialized line drawing instructions. 

Since SEA PLOT was written in as general a form as 
possible, its applications are not limited to the Compass 
Hammer project. With generality in mind, it was prepared 
as a stand-alone program and was not incorporated into one 
of the Compass Hammer computing routines. SEA PLOT, in 
essence, can be utilized by any user desiring a three- 
dimensional attack scenario plot. It need not be constrained 
to the display of attacking aircraft encounters, but could be 
used from the standpoint of a ship defending itself from 
a single aircraft, evaluation of defensive missile design, 
or evaluation of weapon scenario development. Another 
possible use of SEA PLOT would be in conjunction with the 
Pool Input Program (PIP) , which currently only provides 
two-dimensional plotted output. 
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SEA PLOT Output with SCALE= 
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Fij^ure 1-2: SEA PLOT Output with SCALE = 
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i-gure 1-3; SEA PLOT Output with SCALE = 



Detailed discussion of SEA PLOT techniques and use 
follow in Chapters II and II which cover the conic pro- 
jection and program input and organization, respectively. 
Appendix A provides a complete program listing, Appendix B 
contains a listing of the conic projection routines, and 
Appendix C discusses use of the Tektronix developed inter- 
active graphics software PLOT-80 on the Tektronix 4081 
graphics terminal. 
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II. CONIC PROJECTION 



This chapter discusses the projection technique employed 
by SEA PLOT to provide a perspective two-dimensional 
display of an object in three-dimensional space. The conic 
projection technique was selected for SEA PLOT from among 
a multitude of algorithms [Ref. 1]. This technique is 
particularly well suited to the purpose of displaying a 
three-dimensional object as if it were being physically 
viewed by the user's eye in that not only is the solid 
body orientation in space accurately portrayed, but per- 
spective depth and distance from the viewer's eye are 
preserved. Thus, it becomes an excellent tool for viewing 
an aircraft-ship engagement as observed from a distant 
location. 

As with any effort to display a solid object using a 
mathematical definition of body doorcinates, the elimination 
of hidden lines and surfaces as they are cloaked from view 
by surfaces in a more prominent viewing position presents a 
problem to the programmer. Moreover, the computation time 
required to mask hidden surfaces grows with the square of 
the situation complexity [Ref. 2]. However, the use of 
simplistic solid body models as opposed to detailed body 
structures eliminates the need for a comprehensive hidden 
surface removal routine. Since SEA PLOT is intended to 
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provide a symbolic overview rather than a definitive 
display of an attack scenario, this technique has been 
adopted. Consequently, the aircraft body is displayed in 
triangular form approximating the wing planform of a delta 
wing aircraft. The triangle is elongated with an included 
tail section to enhance the perspective orientation of the 
aircraft in three-space. Likewise, the target ship is 
projected as a two-dimensional hull form oriented along 
a user’s input target heading. However, should future use 
of SEA PLOT dictate more comprehensive target and attacker 
models. Reference 3 characterizes ten hidden surface algo- 
rithms which a programmer may employ. 

The conic projection technique involves placement of 
the viewer’s ’’eye” in an earth- fixed coordinate system. 
Associated with the positioning of the eye is a cone of 
vision encompassing the field of view and a projection 
plane onto which the scenario is to be displayed. SEA 
PLOT automatically determines these positions based upon 
the input flight profile of the attacking aircraft. The 
viewing plane is placed 2000 meters beyond (as viewed from 
the origin) the attacker’s initial flight path point, 
oriented perpendicular to a line drawn through the earth- 
fixed origin (the location of the ship) and the initial 
attacker position. The eye coordinates are placed 2000 
meters beyond the viewing plane along the same line. This 
automatic placement of the viewing system eliminates the 
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requirement for the user to determine an optimal viewing 
point which, if ill-placed, could result in plot distortion 
should the field of view exceed the cone of vision. Should 
another viewing position relative to the attacker’s air 
track be desired. Chapter III discusses the procedure to 
alter the perspective position. Appendix B contains the 
actual routines employed by SEA PLOT to accomplish the conic 
projection and can be used as a ’’skeleton program” to be 
applied to three-dimensional graphics endeavors of some 
other nature. 

The conical projection method described in this 
chapter is essentially the same as that used in some of the 
displays with the Differential Maneuvering simulator at the 
NASA Langley Research Center [Ref. 4]. It is summarized in 
Reference 5, from which the following derivations were taken. 
Reference 5 also discusses additional applications of the 
technique . 

First consider an aircraft at some distance from the 
center of an earth-fixed coordinate system. The coordinate 
system, with axes labeled Xe , Ye, Ze is right handed, the 
Xe-Ye plane tangential to the earth at the origin, with the 
Z axis pointing up as shown in Figure 2-1. 

The solid body (in this case the aircraft) center of 
gravity (origin of the aircraft body-fixed coordinates) is 
then positioned at Xe , 7e , Ze , with the body orientation in 
three-space given by the Euler angles ip , d , ip corresponding 
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to heading, pitch, and roll, respectively. The aircraft 
body axes (Xb, Yb , Zb) are chosen such that Xb points from 
the origin toward the nose, with the Yb axis pointing in 
the direction of the left wing which gives a conventional 
right-handed system with the Zb axis pointing up from the 
aircraft center. Coordinates of the aircraft body are 
given by Xp(i), Yp(i), Zp(i) which are relative to the air- 
craft center of gravity. These body coordinates can then be 
transformed to earth-fixed coordinates Xp'(i), Yp*(i), Zp'(i) 
using 
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where |pj is the direction cosine matrix from earth-fixed 
axes to body axes; 
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After transforming the aircraft body coordinates to 
earth-fixed coordinates, the viewing plane must be established 
onto which the conic projection will be made. In so doing, 
it is convenient to establish a new coordinate system 
Xev, Yev, Zev, with the origin at the eye point (Xeye, 

Yeye , Zeye) . The X axis of this coordinate system extends 
from the eye point to some view point located on the viewing 
plane. This, then, represents the axis of the projection 
cone. The viewing coordinate system is obtained by two 
rotations of the earth-fixed system. First, a rotation is 
made about the Ze axis by the angle 9ev, the angle between 
the Xe axis and the projection of the line connecting the 
eyepoint and the viewpoint onto the Xe-Ye plane. The second 
is a pitch rotation about the new Y axis by an angle 9ev. 

This orients the cone of vision in space. From Figure 2-1, 
it can be seen that 



Q "1 ,Zview - Zeye. 

eev = tan ( ftip ) 



where 

R2D = [(Xview - Xeye)^ + (Yview - Yeye)^]^'^^ 



Any point in the earth-fixed coordinate system can 
now be expressed in the viewing coordinate system by 
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Xev 




X'p(i) 




Xeye 


Yev 
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Y'p(i) 
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Zeye 



where [A] is the previously defined [D] matrix with 

0 = 0ev 
\l> jpev 
<j) (j)ev 

The final step, after transforming the body coordinates 
of the solid object to viewing coordinates is to make a 
two-dimensional projection of the three-dimensional object 
onto a flat picture plane. (For SEA PLOT, this picture 
plane is the plotting paper; but more generally, it would 
be the screen of the graphics display console.) To do this, 
the distance from the viewer (user's eye) to the picture 
plane is defined as the distance XPLANE. (See Figure 2-2.) 
This permits determination of the two dimensional plotting 
coordinates [YSCOPE, ZSCOPE) for the projection of a three- 
dimensional object onto a two-dimensional plotting plane. 
From Figure 2-2, it can be seen that the two triangles 
OSF and OP"P* are similar. Therefore, 

XPLANE ^ YSCOPE 
Xev ~ Yev 
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or 



YSCOPE 



Yev 

Xev 



X XPLANE 



Also similar triangles are OPF' and OPP ' . 



ZSCOPE _ OF 
Zev OP^ 



From the previous set of similar triangles 



OF' _ XPLANE 
OP ' Xev 



Hence , 



ZSCOPE 



^ X XPLANE 



Appendix B provides the FORTRAN coding 
these derivations. 



Therefore , 



it follows that 



to implement 
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Figure 2-1: Conical projection of aircraft image 

on to picture plane. 
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Figure 2-2: Transformation from the EV system into the SCOPE system. 
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III. PROGRAM INPUTS AND ORGANIZATION 



This chapter discusses SEA PLOT required input data and 
the organization o£ the program itself. It provides informa- 
tion on the input data and input format as well as detailed 
discussion of the thirteen different routines that comprise 
SEA PLOT. 

A. PROGRAM INPUTS AND FORMAT 

Program inputs for SEA PLOT are intended to provide 
the user with a wide variety of options in designing the 
desired display output. SEA PLOT inputs fall into two 
categories: plotting parameters and flight path coordinate 

points. A maximum of 200 flight path coordinate points may 
be plotted. Following is a description of input elements 
and format. 

1 . Card 1 - Plotting Parameter Input Card 

a. Col 1-3: Enter integer number of flight 

path milestones, right justified. 

b. Col 4: Enter integer value of the type of 

ship defensive missile envelope. The two types of envelope 
display available are a solid cylinder (Figure 3-1) or a 
bold black range ring drawn on the sea surface (Figure 3-2). 
Integer inputs are: 

1 - range ring. 

0 - solid cylinder. 
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c. Col 5: Enter integer value indicating if a 

plot of the coordinate axis is desired as follows: 

1 - axis plot is desired (Figure 3-3) . 

0 - no axis plot desired. 

d. Col 6-15: Enter decimal value corresponding 

to the desired plot scale. Figures 1-1, 1-2, 1-3 show 
representative scaled plots corresponding to SCALE=1 . , 2 . , 3 . , 
respectively. This input parameter essentially magnifies the 
plot with the center of magnification at the origin of 

the earth-fixed coordinate system. No distortion of the 
perspective view results from increased scale factors. 

e. Col 16-25: Enter floating point ship defensive 

missile range in meters. 

f. Col 26-34: Enter floating point value of 

target heading in degrees, counter clockwise from the x axis. 

2 . Cards 2 - 201 (One Card per Milestone) 

a. Col 1-10; Floating point value of aircraft 
X coordinate in meters. 

b. Col 11-20; Floating point value of aircraft 
y coordinate in meters. 

c. Col 21-30: Floating point value of aircraft 

z coordinate in meters . 

d. Col 31-40; Floating point value of aircraft 
roll in degrees. 

e. Col 41-50: Floating point value of aircraft 

pitch in degrees. 
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f. Col 51-60: Floating point value of aircraft 

heading in degrees. 

B. PROGRAM ORGANIZATION 

SEA PLOT consists of thirteen different routines 
which can be categorized as plot management routines, 
graphics routines, conic projection routines, and plotting 
routines. This organization is intended to provide a 
smooth flow of data from one program element to another as 
well as to present a reasonably coherent program structure 
to the user. Sufficient commentary data has been included 
in the program to facilitate future changes or alterations 
in as uncomplicated a manner as possible. Program comments, 
coupled with the detailed discussion in this chapter, should 
provide sufficient program background such that program 
modifications to meet a variety of needs could be relatively 
expeditiously made. However, a word of caution is necessary 
when altering data or variable elements. Much data is 
passed from routine to routine in the form of Common Blocks. 
While this technique provides an efficient data flow, it 
does detract from the comprehensibility of the coding. 
Consequently, altering data in one subroutine could have 
unpredictable consequences if it impacts on calculations 
taking place in other subroutines. The majority of 
the information passed in common involves input data arrays 
and scaling parameters, some from external sources and some 
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internally generated. Before altering these variables, 
one should thoroughly understand their uses in the various 
algorithms in which they are employed. 

1 . Plot Management Routines 

These routines consist of those algorithms required 
to input data, determine plotting scales, label plots, and 
set up the eye-view coordinate system. Three routines com- 
prise this section; they are MAIN, SEAPLT, and AXIS. 

Each of these is described below, 
a. MAIN Program 

MAIN is straightforward in that its primary 
purpose is to read the input data in the format previously 
described. It also establishes the common blocks "AIR," 
"PARAM,'* and "INPUT.” Variables and their descriptions 
are as follows : 

(1) XAIR, YAIR, ZAIR - Floating point arrays 
of 200 elements each which contain the attacking aircraft 
X, y, z coordinates, respectively. 

(2) ACROLL, ACPTCH, ACHEAD - Floating point 
arrays of 200 elements each which contain the aircraft roll, 
pitch, and heading corresponding to each coordinate position. 

C3) MDISP - Integer input value determining 
the type of missile envelope display to be plotted. 

MDISP=1 causes a single bold range ring to be plotted 
on the surface (z=0) . MDISP=0 causes a three-dimensional 
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cylinder to be plotted. The height of this cylinder 
corresponds to twice the maximum aircraft input altitude. 
The radius of the envelope is the missile slant range 
which is an input variable. (See Figures 3-1 and 3-2.) 

(4) lAXIS - An integer input value used to 
determine if a plot of the coordinate axes is desired. 
IAXIS=0 indicates no axes are to be drawn. IAXIS=1 
causes the coordinate axes x, y, z to be displayed. 

(5) SCALE - An input floating point value 
used to determine the magnification factor to be applied 
to the plotted output. This permits the user to "zoom in" 
on the center of the plot, magnifying it by the value of 
SCALE. Figures 1-1, 1-2 and 1-3 demonstrate the use of 
this parameter. 

(6) RNGMIS - An input floating point variable 
used to define the maximum range of the ship’s defensive 
missiles . 

(7) THEAD - An input floating point variable 
corresponding to the heading of the ship target measured 
in degrees from the axis. 

Printed output from MAIN is a list of input 
variable values and aircraft flight path data, 
b . SEAPLT 

SEAPLT is the driving routine for the entire 
program. It searches for maximum flight path values, 
positions the eye and viewplane coordinates, and sequences 



28 






I; 



I: 



\ 





through the plotting routines to cause the graphic output. 

The viewing plane is positioned 2000 meters behind the 
initial flight path position along a three-dimensional line 
from the origin through the initial position. The center 
of the eye-view coordinate system is placed on this same 
line 2000 meters beyond the viewing plane. Because of the 
low aircraft altitudes employed in current war-at-sea tactics, 
this routine performs an altitude enhancement which adjusts 
aircraft input altitudes such that the maximum plotted 
altitude corresponds to the viewing height. This preserves 
a perspective in which the viewer is always looking along 
or down upon the flight path. In that all of the flight 
path altitudes are increased by this same value, no plot 
distortion occurs, overall perspective is maintained, and 
altitude maneuvers are enhanced. 

Variables used in SEAPLOT are as follows: 

(1) XVIEW,YVIEW,ZVIEW - x, y, z coordinates of 
the center of the viewing plane. 

(2) XEYE , YEYE , ZEYE - x, y, z coordinates (earth- 
fixed) of the eye position. This becomes the origin of 

the eye-view coordinate system. 

(3) XPLANE - Physical distance from the viewer 
to the plotting surface. This variable is particularly 
germaine to plotting on a graphics terminal. Here it is 
fixed at 30 cm. 



29 



(4) DX - Separation of range rings from the plot 
origin. DX is currently fixed at 5000 meters. 

(5) BESTZ - Altitude enhancement factor applied 
to the input aircraft flight profile altitude. 

(6) XMAX - Maximum aircraft x coordinate. 

(7) YMAX - Maximum aircraft y coordinate. 

(8) ZMAX - Maximum aircraft z coordinate. 

(9) THETA - Angle in radians of the line drawn 
through the origin and the first aircraft coordinate position. 

c. AXIS 

In addition to drawing coordinate axes, subroutine 
AXIS also performs plot labeling and display scaling. Con- 
sequently, this subroutine must be called whether or not a 
plot of the coordinate axes is desired. The only parameter 
required by AXIS is NFLAG. This is the lAXlS argument pre- 
viously described. 

Scaling performed by AXIS is accomplished in such a 
manner as to center the plot on a Versatec 20-inch by 20-inch 
plotting sheet. To do this, AXIS first projects a circle of 
radius corresponding to the length of a vector drawn from 
the earth-fixed origin out to the point given by the coordi- 
nates ( XMAX, YMAX ) . This circle consists of 180 data points. 
After calling the conic projection routines to make the 
projection circle, AXIS searches for the longest vector from 
the origin out to the periphery of the circle. Since the 
projected circle will appear as an ellipse whose elongation 
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will depend on the viewing angle, this longest vector will 
vary from view to view. Once this longest projection vector 
is known, AXIS reorigins the plot to the center of the 
plotting sheet at plotter coordinates (10,10), then scales 
the length of the longest projection vector, named SFACT, 
such that its physical size is 10 inches on the plotter. 

This permits the development of an over-all plot scaling 
factor XSCALE where: 



XSCALE = 10.0/ SFACT 

Plot magnification can be accomplished to meet the 
user's specifications by multiplying XSCALE by the input 
variable SCALE. This, of course, will result in the ultimate 
computation of many vectors whose length will exceed the 
dimensions of the plotting paper. However, these vectors 
will be clipped at the edge of the plotting sheet and will 
not degrade the plotted output; though a diagnostic, which 
may be ignored, will be generated on the ouptut listing 
containing the flight path data. 

The final function of AXIS is to label the plot. 
The labeling instructions in the program are software 
dependent. Those in AXIS are written specifically for 
VERSATEC software. Hence, careful attention must be paid 
to this routine when adapting it to another graphics system. 

Axis options and labeling are shown in Figure 3-3 
which shows the axis output. 
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2 . Graphics Routines 



The graphics routines described in this section are 
those elements which accomplish the actual pen moves required 
to draw lines or dots on the plotting sheet. Thus, these 
routines are highly software dependent and represent the 
major changes required to adapt SEA PLOT to another graphics 
system. SEA PLOT requires only three types of graphics 
operations: drawing solid lines, drawing vertical dashed 

lines, and plotting dots. The two subroutine descriptions 
which follow discuss these operations, 
a. Subroutine DRAW 

This subroutine causes an array of projected 
coordinate data points to be connected either by a solid 
line or to be plotted as dots on the plotting sheet. Up 
to 200 data points may be plotted. Projected coordinates 
are first scaled to a 20-inch by 20-inch plot size using the 
XSCALE factor developed in AXIS. They are then shifted 
from the earth-fixed coordinate center to the plotting sheet 
coordinate center by adding a value XORIG, or YORIG which 
corresponds to the displacement of the projected earth 
origin from the plotter origin. DRAW next causes a pen up 
move to be made to the initial array position and then either 
draws a line connecting the points or plots a dot correspond- 
ing to each coordinate position. 
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Variables used in DRAW are: 



(1) X,Y - 200 element input arrays of projected 
data to be plotted. 

(2) XPLOT, YPLOT - 200 element arrays containing 
the scaled projection data to be plotted. 

(3) N - Integer value indicating the number of 
coordinate pairs to be plotted. 

(4) ITYPE - Integer value which determines the 
type of plot to be drawn. ITYPE=0 causes dots to be plotted. 
ITYPE=1 causes solid lines to be drawn. 

b. Subroutine VDASH 

Subroutine VDASH causes a vertical dashed line, 
parallel to the z axis, to be drawn between two points whose 
coordinate values are passed as arguments. It should be 
emphasized that only a vertical line is drawn. Should the 
two points to be connected by the line not be vertically 
positioned, the routine will draw a dashed line from the 
first point up to the Y value of the second point. As with 
subroutine DRAW, all data values are scaled to the XSCALE 
factor generated in AXIS. 

Variables used in VDASH are: 

(1) XI, Y1 - Projection coordinates of the lower 

data point. 

(2) X2,Y2 - Projection coordinates of the upper 

data point. 



points . 



(3) DY - vertical separation between the data 
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3 . Conic Projection Routines 



The three subroutines that follow are those routines 
used to implement the conic projection technique derived in 
Chapter II. They are also reproduced in Appendix B to 
serve as a ready reference "skeleton program" to be used in 
other graphics efforts. 

a. Subroutine MULT 

MULT performs a matrix multiplication of a 
square 3x3 matrix with a column vector, returning the 
resultant column vector to the calling program. The opera- 
tion performed is 







— — 


Xnew 




X 


Ynew 


= [D] 
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Znew 




Z 


— ^ 




^ 



where the vector [X,Y,Z] and the matrix [D] are passed to 
the subroutine as arguments. 

b . Subroutine DIRCOS 

DIRCOS is a rather simplistic routine which 
merely sets up the direction cosine matrices discussed in 
Chapter II. Inputs to the routine are the Euler angles 
psi, theta, and phi corresponding to solid body roll, pitch, 
and yaw in three-space. DIRCOS output is the three-dimensional 
matrix jpj or employed by the projection. Also available 
as an output is the transpose of [d1. In that the coding used 
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in DIRCOS is a straightforward assignment of variables to 
matrix positions, no further discussion will be attempted 
here . 

c. Subroutine CONPRO 

CONPRO is the primary routine for the conic 
projection scheme implementing the derivations in Chapter II. 
Inputs to CONPRO are the center of gravity coordinates of 
the solid body to be projected, the number of defined body 
coordinates, the array of body coordinate points, and the 
body Euler angles in degrees measured counterclockwise from 
the earth-fixed X axis. After performing all translations 
and rotations required by Chapter II, CONPRO outputs the 
plotting surface coordinates of the projected solid body. 
These coordinates are passed by the calling program to the 
DRAW routine where they are scaled and plotted. 

Variables used in CONPRO are: 

(1) CG - Three element array containing the 
X,Y,Z earth-fixed coordinates of the solid body center of 
gravity . 

(2) NPOINT - Integer value indicating the 
number of data points used to define the solid body. 

(3j XP - Array of body X coordinate values in the 
body-fixed reference system. This array is variably dimen- 
sioned to NPOINT elements. 

(4) YP - Array of body Y coordinate values in the 
body-fixed reference system. This array contains NPOINT 
elements . 
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(5) ZP - Array o£ NPOINT body Z coordinate values 
in the body- fixed reference system. 

(6) EULER - Three element array containing the 
Euler angles for roll, pitch, and heading in degrees. 

(7) YSCOPE - Output array of NPOINT projected 
plotter X coordinates. 

(8) ZSCOPE - Output array of NPOINT projected 
plotter Y coordinates. 

4 . Plotting Routines 

The five subroutines that follow are those algorithms 
which develop the different picture segments that comprise 
the SEA PLOT output. The five elements displayed are: the 

circular surface grid consisting of 5000-meter range rings, 
the aircraft ground and air tracks, the aircraft body plot, 
the target body plot, and the defensive missile launch 
envelope. Since more sophisticated graphics languages such 
as PLOT 80 employ the concept of graphics control blocks 
(discussed in Appendix C) or picture segments, SEA PLOT 
was organized along these same lines to facilitate easier 
transition to one of these more powerful software systems. 
Using this scheme, each subroutine will cause a specific 
picture segment to be displayed as directed by the sequencing 
routine SEAPLT. As a result, the user has the option of 
deleting those picture elements which are not desired. For 
example, by circumventing the call to the target subroutine, 
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one could get a plot of only the surface grid and aircraft 
flight profile without a target being shown, 
a. Subroutine CGRID 

Subroutine CGRID projects and draws a series of 
concentric range circle centered at the target position at the 
origin. As presently written, these range circles are drawn 
at 5000-meter increments as indicated by the variable DX 
discussed in SEAPLT. Each range circle is defined by plane 
polar coordinates consisting of 181 data points per circle. 

The number of circles drawn is a function of the maximum 
X coordinate for the aircraft flight profile such that the 
number of circles is the quotient of the maximum X range 
and the radius increment DX. 

Variables used in CGRID are: 

(1) X,Y,Z - 181 element arrays containing the 
unprojected coordinates of the data points defining each 
range ring. The Z array is set to zero corresponding to 
zero altitude. 

(2) YSCOPE , ZSCOPE - 181 element arrays containing 
the projected plotter coordinates of the range circle data 
points. These arrays are passed to the DRAW subroutine where 
they are scaled and plotted. 

(3) DX - Interval between range circles in the 
earth-fixed reference system. 

(4) RADIUS - Radius of each range circle in the 
earth-fixed system. 
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Figure 3-4 portrays the output from CGRID. 



b. Subroutine ACPLOT 

Subroutine ACPLOT plots up to 200 aircraft 
flight path coordinate points in a three-dimensional 
perspective view. To enhance the view, both ground and air 
tracks are displayed. The ground track is plotted by setting 
the aircraft Z coordinate to zero, then projecting the X,Y,Z 
coordinates and plotting the data as series of dots corre- 
sponding to each input point. In a like manner, the air 
track is projected without zeroing the Z coordinate. A 
solid, bold line connects all air track data points. To 
reinforce the altitude perspective, an optional dashed 
vertical line is drawn between corresponding air and ground 
track data points. This vertical line option parameter is 
a program change option which is specified by the calling 
program not by a user input. 

Aircraft maneuvers in roll, pitch, and yaw are 
displayed by plotting the aircraft body at positions corre- 
sponding to every 2000 meters of ground travel in the earth- 
fixed system. This is accomplished by a check of ground 
distance traveled between input points followed by a call 
to subroutine ACRAFT after 2000 meters of displacement has 
been detected. Because the surface area of the aircraft 
body is insignificant with respect to the total surface area 
displayed by the plot, a magnification factor is applied 
to the aircraft body. It has been experimentally determined 
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that a magnification factor of 15 permits an optimal body 
display. 

Inputs to ACPLOT include the number of milestones 
to be plotted, the vertical line option parameter, and 
magnification factor. These and the other internal variables 
used in this routine are discussed as follows: 

(1) N - Integer number of aircraft milestones 
to be plotted. 

(2) VLINE - Decimal flag used to exercise the 
dashed vertical line option. If VLINE=1, the vertical 
altitude lines are drawn. Otherwise, no vertical lines 
are drawn between tracks. 

(3) AMAG - Decimal magnification factor to be 
applied to the aircraft body plot. 

(4) X,Y,Z - 200 element arrays containing the 
earth-fixed center of gravity coordinates (input milestones) 
for the body plot and the flight profile tracks. 

(5) XPLOTA,YPLOTA - 200 element arrays con- 
taining the projected X,Y coordinates for the air track. 

(6) XPLOTG ,YPL0TG - 200 element arrays con- 
taining the projected X,Y plotter coordinates for the ground 
track. 

Figures 3-5 and 3-6 show the plot output options 

for ACPLOT. 
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c. Subroutine ACRAFT 

ACRAFT contains the predefined aircraft body 
coordinates in the body-fixed coordinate frame. This routine 
applies the magnification factor to the body coordinates, 
causes a projection of the body to be made and a plot output 
to be generated. The actual body coordinates are defined 
in data statements in the program and are not user verifiable. 

Inputs to ACRAFT are the Euler of angles roll, 
pitch and yaw, the center of gravity coordinates in the earth- 
fixed reference, and the magnification factor. These and 
other variables used by ACRAFT are defined as follows: 

(1) EULER - A three element array containing 
the body roll, pitch, and yaw angles in degrees. 

(2) CG - A three element array containing the 
X,Y,Z earth-fixed coordinates of the center of gravity. 

(3) XAC, YAC, ZAC - Eight element arrays 
containing the predefined aircraft body coordinates in the 
body reference frame. 

(4) AMAG - Magnification factor to be applied 
to the body coordinates. 

The output from ACRAFT can be seen in Figures 

3-5 and 3 - 6 . 

d. Subroutine TARGET 

TARGET performs essentially the same function as 
ACRAFT except that it contains the predefined body coordinates 
of the ship target. Because an excessive number of plotted 
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vectors tends to make the overall plot output too busy, the 
target is displayed as a hull planform only. The planform 
is shaded to accentuate its position on the plot. At present, 
the target is plotted at the earth-fixed origin. However, 
this is not a rigid requirement since this particular 
routine has been written in a general manner to permit 
positioning the target anywhere in the viewing area. Thus, 
the location of the target center of gravity is an input 
from SEAPLT, the calling program, and can be altered to fit 
the user's needs. The other input to TARGET is the ship 
heading measured in degrees from the earth- fixed X axis. 

This also is a user specified input parameter. 

Variables used in TARGET are; 

(1) CG - A three element array containing the 
target earth-fixed center of gravity coordinates. 

(2) THEAD - Target heading measured in degrees 
counterclockwise from the earth- fixed X axis. 

(3) TMAG - Magnification factor applied to the 
body fixed target coordinates. This variable is set to the 
same value as AMAG in ACRAFT. 

(4) XTGT ,YTGT, ZTGT - 26 element arrays containing 
the predefined body coordinates of the target in the body- 
fixed reference system. 

(5) XSCOPE , YSCOPE - 26 element arrays containing 
the projected coordinates of the target in the plotter refer- 
ence system. These are scaled and plotted by a call to DRAW. 
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Figures 3-7 and 3-8 demonstrate the plot 
options available to TARGET. 

e. Subroutine LAUNCH 

LAUNCH causes a plot display of the target 
ship's defensive missile envelope to be drawn. Two options 
are available. They are a three-dimensional cylinder or a 
bold black circle plotted on the surface as shown in Figures 
3-1 and 3-2, respectively. Selection of the desired 
envelope plot is a user input parameter. Inputs to LAUNCH 
are the missile range and the display type. 

The cylindrical display is defined using 
cylindrical coordinates such that the height of the 
cylinder cooresponds to twice the maximum aircraft input 
altitude. Cylinder construction consists of ten circles 
vertically separated by a distance of one fifth the maximum 
aircraft altitude. Each circle consists of 180 dots. 

The range ring display merely consists of 180 
coordinate points expressed in plane polar coordinates. The 
radius of the circle corresponds to the maximum missile 
launch range. 

Variables used in LAUNCH are: 

(1) RANGE - Maximum missile launch range. This 
is an input parameter from the calling program (SEAPLT) . 

(2) IDISPL - Integer display option flag such 
that IDISPL=1 causes the single range ring display. Other- 
wise the cylinder display is drawn. 
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rigui'e 3-1; Cylinier I.auijch Envelope Display 
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Figure 3-2: R-^ngf- Rir.g Liunch Fnvelope Display 
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Figure 3-3: Coordinate Axis Display Option 

Output from Subroutine AXIS 
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Figure 3-4; CGKIl) Range Ring Output 
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Figure 3-5: ACPLOT Output without Altitude 

Enhancement Lilies 
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ri^ur-'r 3-6: A/'FLOT Output v;ith Altitude 

Er.lianceinent Lines 
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Figur-e 3-7: TARGET output. Ship heading = 000 
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IV. SUMMARY AND CONCLUSIONS 



SEA PLOT is a three-dimensional computer graphics 
program written to display a war-at-sea encounter between a 
single ship and a single aircraft, in which the ship is 
assumed to be the target and is located at the origin of an 
earth-fixed coordinate system. Roll, pitch, and yaw maneuvers 
of the attacking aircraft are shown along a three-dimensional 
flight path which depicts the attacker's coordinate positions 
in the earth-fixed reference system as well as the attitude 
of the aircraft body in three space. Up to 200 flight path 
coordinates can be plotted. To accomplish the projection of 
a three-dimensional scenario onto a two-dimensional plotting 
surface, a conic projection technique is used in which a 
viewing plane and an eye position are automatically computed 
based upon the attacking aircraft's initial position. 

Displays available to the user include: (1) the earth- 

fixed coordinate axes, (2) 5000-meter circular range rings 
centered on the target, (3) ship defensive missile launch 
envelope displayed in solid cylindrical or plane circular 
form, (4) aircraft air track, (5) aircraft ground track, 

(6) aircraft maneuvers of roll, pitch, and yaw, (7) ship 
target body display located at the earth-fixed origin. To 
generate these displays, SEA PLOT was written using thirteen 
subroutines which fall into four functional categories: plot 
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«IV 



management routines, graphics (vector drawing) routines, 
conic projection routines, and plotting (picture segment) 
routines. This organization was selected to allow for easy 
conversion from the present VERSATEC graphics software to 
other graphics languages such as the Tektronix developed 
PLOT 80. 
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SEA PLOT PROGRAM LISTING 
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READ (5,900) NPOINT ,MDISP , lAXIS , SCALE ,RNGMIS ,THEAD 
PRINT FLIGHT PROFILE PARAMETERS 
WRITE (6,903) 

WRITE (6,904) SCALE, RNGMIS ,THEAD, NPOINT 



WRITE (6,905) 

READ AND PRINT FLIGHT PROFILE 
DO 10 I=l,NPOINT 



to 
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POSITION EYE POINT 2000M BEYOND VIEW POINT - THIS 
ORIENTS THE VIEWING DIRECTION 
EYERAD=VRAD+2000 . 

XEYE=EYERAD*COS (THETA) 

YEYE=EYERAD*SIN (THETA) 
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CALL LAUNCH (RNGMIS ,MDISP) 
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CALL PLOT (XABS(3) ,YABS(3) ,+2) 
CALL WHERE (YLABX , YLABY , DFACT) 
CALL PLOT(0.0,0.0,+3) 

CALL PLOT(XABS(4) ,YABS(4) ,+2) 
CALL WHERE (ZLABX , ZLABY , DFACT) 
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RETURN 

END 



SUBROUTINE DRAW (X , Y ,N , ITYPE) 

DIMENSION X(N) ,Y(N) ,XPLOT(200) ,YPLOT(200) 

COMMON/ PSCALE/ XVAL , YVAL , XV , YV , XSCALE , YSCALE 
COMMON/ GRI D/ XMI N , XMAX , YMI N , YMAX , DX , ZMAX , XORI G , YORI G 
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DRAW PLOT USING LINES 
30 DO 40 1=2, N 

CALL PL0T(XPL0T(I) ,YPL0T(I3 ,+2) 
CONTINUE 
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I 



XNEW=D(1,1)*X +D(1,2)*Y +D(1,3)*Z 
YNEW=D(2,1)*X +D(2,2)*Y +D(2,3)*Z 
ZNEW=D(3,1)*X +D(3,2)*Y +D(3,3)*Z 
RETURN 
END 
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SUBROUTINE DlRCOS(PSI , THETA, PHI ,D,TRANSP) 
DIMENSION D(3,3) 

SUBROUTINE COMPUTES DIRECTION COSINE MATRIX D 
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SUBROUTINE CONPRO (CG , XP , YP , ZP ,NPOINT , EULER, YSCOPE , ZSCOPE) 
DIMENSION CG(3) ,XP(NPOINT) , YP (NPOINT) , ZP (NPOINT) ,EULER(3) 
DIMENSION YSCOPE (NPOINT) , ZSCOPE (NPOINT) ,A(3,3) ,C(3,3) 
COMMON/EVIEW/XVIEW , YVIEW, ZVIEW,XEYE , YEYE , ZEYE , XPLANE 
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SUBROUTINE MAKES A CONIC PROJECTION OF POINTS IN A 3 SPACE PASSED 
IN XP,YP,ZP. PROJECTION BASED ON VIEWING COORDINATES, CENTER OF 
GRAVITY, AND EULER ANGLES (ROLL , PITCH , HEADING) . XPLANE IS THE 
DISTANCE FROM VIEWER TO THE PROJECTION PLANE. 
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TRANSLATE TO ORIGIN OF EYE -VIEW SYSTEM 
XPR=XPR+CG(1) -XEYE 
YPR=YPR+CG(2) -YEYE 
ZPR=ZPR+CG(3) -ZEYE 
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Z(I)=0.0 

CONTINUE 
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SUBROUTINE ACPLOT (N , VLINE , AMAG) 

COMMON/ AI R/XAI R , YAI R, ZAI R , ACROLL , ACPTCH , ACHEAD 
DIMENSION XAIR(200) ,YAIR(200) ,ZAIR(200) 
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CALL NEWPEN(S) 

CALL DRAW(XPLOTG,YPLOTG,N,0) 



DEFINE MAGNIFICATION VALUE FOR DISPLAY OF A/C BODY 
CALL NEWPENC2) 

SET UP MILESTONE COUNTER FOR DISPLAY OF A/C BODY 
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PREDEFINED AIRCRAFT BODY COORDINATES: 
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SUBROUTINE TARGET PROJECTS AND PLOTS A GENERIC SHIP IMAGE CENTERED 
AT COORDINATES PASSED IN CG ARRAY, WITH HEADING PASSED IN THEAD. 
IMAGE IS MAGNIFIED BY A FACTOR OF TMAG . TARGET BODY COORDINATES ARE 
PREDEFINED IN THIS ROUTINE. 



AS WRITTEN, THIS ROUTINE DRAWS ONLY THE LOWER PLANE OF THE HULL. 
THIS IS TO RELIEVE A CROWDED PLOT. HOWEVER, IF A HIDDEN LINE 
ROUTINE IS INCLUDED WHICH WILL CLEAN UP THE PLOT, THE REST 
OF THE BODY COORDINATES ARE PROVIDED IN THE DATA STATEMENTS. 
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AND ORIGIN 

DO 20 1=1,6 

XSC0PE(I)=XSC0PE(I)*XSCALE+X0RIG 
YSCOPE (I)=YSC0PE (I) *YSCALE+Y0RIG 
20 CONTINUE 



CALL TONE(0. ,0. ,IPAT,-5) 
CALL TONE (XSCOPE , YSCOPE , 
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CALL DRAW(XPLOT,YPLOT, 181,1) 

80 RETURN 
END 
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THETEV=ATAN4 (YARC ,R2D) 
CALL DIRCOS(PSIEV,THETEV 
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CPSI=COS(PSI) 
STHET= -SIN (THETA) 
CTHET=COS (THETA) 
SPHI=SIN(PHI) 
CPHI=COS(PHI) 
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I. INTRODUCTION 



This appendix discusses the Tektronix graphics software 
package PLOT 80 and its use on the Tektronix 4081 graphics 
computer. Because Reference 6 provides an adequate dis- 
cussion of the Tektronix 4081 equipment and system operation, 
computer light-off and operation will not be reiterated 
here. Instead, discussion will be focused on the graphics 
language PLOT 80 and its uses in the Tektronix 4081 program 
mode . 

PLOT 80 is a particularly powerful graphics tool con- 
sisting of 374 different subroutines covering a wide range 
of graphics concepts including bit stream manipulation, 
vector drawing, test input and output, and picture display. 
These subroutines are discussed in detail in the PLOT 80 
user’s manual. Reference 7. Because PLOT 80 represents a 
complex, sophisticated software system. Reference 7 should 
be reviewed by a prospective user prior to any attempt at 
programming of the Tektronix 4081. This appendix will 
discuss some of the more basic concepts employed by PLOT 80 
and provide the programmer some useful background informa- 
tion which is intended to complement the information in 
Reference 7. 

PLOT 80 will be discussed here in the context of the 
4081 program mode of operation as opposed to the host mode 
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discussed in Reference 6. In the program mode, all computa- 
tions are accomplished using the hardware organic to the 
4081 system itself; unlike the host mode in which a host 
computer, such as the IBM 3033 or PDP-11, performs all calcu 
lations and uses the 4081 for graphics display only. The 
word size of the 4081 does inhibit many mathematical computa 
tions , such as the transforms used in SEA PLOT, which result 
in underflows and hence limits to some extent the program 
mode applications of the 4081. However, this mode does 
provide an excellent tool for investigating graphics 
structures such as the target and aircraft body representa- 
tions used in SEA PLOT. Section II of this appendix dis- 
cusses some of the important basic concepts of PLOT 80. 
Section III will describe the procedures for file creation, 
compiling, linking, and execution on the Tektronix 4081. 
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II. PLOT 80 



The 374 distributed graphics support subroutines (DGSS) 
in PLOT 80 are FORTRAN callable routines that permit the 
user to: 

(1) Perform sophisticated graphics in integer or 
floating point format. For example, routines that perform 
graphic transformations or create polygons, circles, and 
splines are included in the software package. 

(2) Perform graphic input (GIN) in windowspace 
or viewspace coordinates . 

(3) Perform file manipulation. 

(4) Communicate with the host computer. 

(5) Perform assembly language operations including 
bit and character string manipulation, primitive I/O, 

and dynamic memory allocation. 

Before performing any DGSS graphics, a program must 
include the two following steps: 

CALL BFDEV (logical unit, device). 

CALL GONEW (logical unit) . 

The first routine assigns a logical unit (to be defined 
in following paragraphs) to a graphic output device; i.e., 
the display controller. The second routine creates a data 
structure called the graphics control block (GCB) for the 
specified logical unit. By means of the logical unit. 
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therefore, a GCB is linked with a graphic output device as 
shown in Figure A-1. 




Figure A-1 



For example, if the programmer decides that the display- 
controller is to be assigned to logical unit 0, he could 
create a GCB as follows: 

CALL BFDEV (0, ’DC: ') 

CALL GONEW (0) 

He is now ready to create the graphics control block for 
his particular graphics needs. 

The graphic control block contains a sequence of 
parameters that determine exactly where on the surface of 
the graphic output device the graphics information defined 
in a program is to be displayed. The GCB is essential to 
all DGSS graphics. Attempting to output graphics information 
to a device without first creating a GCB causes an error 
condition to be signaled which will terminate execution. 
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No DGSS graphics programiner can entirely ignore the 
concept of a GCB. The GONEW routine, however, allows a 
programmer to get directly to the problem of the graphics 
application without spending much time and space on 
creating a GCB and defining each and every one of its 
parameters . 

The GONEW routine automatically allocates the memory 
required for a GCB and initializes the contents to a set of 
default values. The programmer does not need to know where 
in memory the GCB is located or how much space it occupies. 
Table 14-1 of Reference 7 lists the GCB parameters and the 
default settings. 

Data transfers under the graphics operating system are 
always performed through a logical unit assigned to a 
device. The programmer communicates with the logical unit 
rather than the device itself. This process greatly simpli- 
fies the I/O operation. The programmer does not need to 
do "housekeeping” tasks, such as addressing the device and 
checking the status of the device to see if it is busy or 
ready to send and receive data. All housekeeping is per- 
formed by the proper device driver which is called into 
action when a device is assigned to a logical unit. There 
are 16 logical units, numbered 0 to 15. Initially, all units 
are assigned to the keyboard unit until a BFDEV call is made. 
Any logical unit may be assigned to any device. Graphic out- 
put devices are : 
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(1) display controller 

(2) plotter 

(3) a file 

Another important concept available to the programmer 
is that of a picture segment. Using this concept, the over- 
all plotted output is divided into collections of vectors 
called picture segments. Each segment is numbered and can 
be specified independently of the other segments. The pro- 
grammer can therefore manipulate each segment independently 
of the others. (Recall that SEA PLOT was organized along 
these lines with the coordinate axes, surface grid, launch 
envelope, and flight path each being a separate picture 
segment capable of individual manipulation.) 

Picture segments are also sometimes called refresh 
objects, which refers to how the segment is displayed on the 
screen. Segments can be displayed in either refresh or 
storage. IVhen a segment is drawn in storage, the segment 
is drawn once on the screen and remains visible until the 
screen is erased. A segment displayed in refresh, however, 
is not "stored” on the screen. The segment must be con- 
tinually drawn and redrawn on the screen to remain visible 
(a function automatically carried out and not a programmer 
responsibility) . 

Refresh graphics allow the modification and movement 
of segments without erasing the screen. The effect is 
similar to animation in that a picture is drawn at one point 
on the screen, then redrawn again at another, then again at 



82 



another, and so on. The picture seems to be moving across 
the screen. 

There are three steps in creating a picture segment: 

(1) Initialize the segment assigning it a number 
(up to 255) . 

(2) Specify vector attributes (coordinates, intensity, 

etc. ) . 

(3) Terminate segment creation. 

These are accomplished by the OPEN and CLOSE subroutines. 

To illustrate the use of the GCB and picture segments, 
a possible animation routine for SEA PLOT would be as follows: 
CALL BFDEV (0, 'DC: ') 

CALL GONEW(O) 

C. DRAW COORDINATE AXES (PICTURE SEGMENT 1) 

CALL OPEN (1) 

CALL AXIS (1) 

CALL CLOSE (1) 

C. PLOT AXES IN STORAGE 
CALL FIX (1) 

C. DRAW RANGE RINGS (PICTURE SEGMENT 2) 

CALL OPEN (2) 

CALL CGRID 
CALL CLOSE (2) 

C. PLOT RANGE RINGS IN STORAGE 
CALL FIX (2) 
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C. DRAW AIRCRAFT (SEGMENT 3) IN REFRESH MOVING FROM 
C. MILESTONE TO MILESTONE 
DO 10 I=l,NPOINT 
READ X,Y,Z, ROLL, PITCH, HEAD 
CALL OPEN (3) 

CALL ACRAFT ( Z , Y , Z , ROLL , PITCH , HEAD) 

CALL CLOSE (3) 

C. PLOT SEGMENT 3 IN REFRESH 
CALL POST (3) 

C. TIME DELAY CORRESPONDING TO AIRCRAFT VELOCITY 
CALL TDELAY 

C. ERASE PREVIOUS BODY PLOT 
CALL UNPOST (3) 

C. ITERATE THROUGH REMAINING MILESTONES 
10 CONTINUE 
STOP 
END 
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III. FILE CREATION AND EXECUTION 



This final section discusses Tektronix 4081 program mode 
procedures to create, compile, link-edit, and execute a 
program. Unlike the job control language options of the 
NFS IBM 3033 main frame computer, compiling, linking, and 
executing jobs on the 4081 are entirely distinct operations 
which must be individually performed by the operator. 

A. FILE CREATION 

Creating a local file on the 4081 is accomplished by 
the command "EDIT filename. file type." This operation is 
essentially the same, using the same functions, as the EDIT 
system used on the main frame. When creating a multi- 
structured file, one must first create a library file. 

This library file will then contain each of the individual 
subroutines comprising the overall program. Moreover, each 
subroutine must be individually compiled. 

Using SEA PLOT as an example, one would create a 
library file named SEAPLT containing all of the applicable 
subroutines as follows: 

FORMAT SEAPLT. LIB (4) 

This creates a SEAPLT library with four data blocks 
available for use. 

To create the individual subroutine files in the 
SEAPLT library, the programmer must edit, or create, a 
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separate file in the library for each subroutine as indi- 
cated below: 

EDIT SEAPLT/MAIN.FOR 
INPUT 

enter main program fortran code 

EDIT SEAPLT/CONPRO.FOR 
INPUT 

enter CONPRO fortran code 

EDIT SEAPLT/TARGET.FOR 
INPUT 

enter TARGET fortran code 

or so on until all subroutine files have been entered. 

B. COMPILING 

Each file structure in the library file must be 
individually compiled in the form: 

FORT library name/file name. OBJ, 

LST=library name/f ilename . FOR 

Compiler errors will be in the LST file, which can be 
read by typing the command 

PRINT LST. LST 

Because this becomes somewhat tedious for a file 
with several structures , one can create a batch routine 
which will compile the various subroutines. The following 
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coding will create a compile routine for a batch processor 
j ob . 

EDIT COM. BAT 
INPUT 

FORT LIBRARY NAME/<1> . OBJ , 

LST=LIBRARY NAME/<1>.F0R 
EXIT 

Using this routine, successive subroutines may be rather 
more expeditiously compiled by submitting the batch job 
BATCH COM. filename 

which will compile a single subroutine and must be repeated 
for each subroutine in the library. 

C. LINK-EDITING 

After a program is compiled, the program must be link- 
edited by typing the following commands: 

LINK 

OUT library name/ RUN. OBJ 
LINK library name/main program name 
EDIT library name 
EDIT SYS:DGL,SYS:RTL 

Should the output of the link editor indicate that there 
are undefined terms in the program, this process must be 
repeated until all undefined terms have been resolved. This 
sometimes involves repeating the link-edit step four or 
five times. 
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D. EXECUTION 



Executing the successfully compiled and linked program 
is, refreshingly, a very straightforward procedure. One 
merely enters the command 

USR: library name/RUN 

If any execution errors are detected, they will be dis- 
played on the screen along with any graphics output which 
has been generated before the error was detected. 
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